Ubiquitous V-shape density of states in mixed state of clean limit type II 
superconductors; STM experiment and theory 
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It is demonstrated theoretically and experimentally that the low energy density of states N(E) is 
described by a singular V-shape form N(E) = Nq(H) +a\E\ + 0(E 2 ) for all clean superconductors 
in a vortex state, irrespective of the underlying gap structure. The linear term a\E\ which has 
been not recognized so far is obtained by exactly evaluating the vortex contribution. Based on 
microscopic Eilenberger theory N(E) is evaluated for isotropic gap, line and point node gaps to 
yield a V-shape N(E). STS-STM experiments on NbSe2 and YM2B2C give direct evidence for this. 
We provide arguments on significance of this finding and on relevance to other experiments. 

PACS numbers: 74.25.Op, 74.25. Bt, 74.25.Jb 

Much attention has been focused on exotic super- 
conductors, ranging from high T c cuprates, Ce and U 
based heavy Fermion materials, filled skutterdites such 
as PrOs 4 Sbi2 to cobaltites Na^CoOryHaOQ, Q. The 
identification of the Cooper-pair symmetry consists of 
two parts; its parity and the gap structure. The former 
is responsible for the spin structure of a pair which is 
either singlet or triplet. This can be probed by directly 
measuring the spin susceptibility through NMR Knight 
shift experiment under an applied field. The latter gap 
structure is related to the orbital symmetry of a Cooper 
pair. This can be probed by thermodynamic measure- 
ments via a variety of experimental methods, such as 
temperature (T) dependence of specific heat C(T), ther- 
mal conductivity k(T) or nuclear relaxation time T\(T) 
in NMR-NQR experiments. 

The basic principle of this identification for the gap 
structure is based on the fact that the energy (E) de- 
pendence of the density of states (DOS) N{E) near the 
Fermi level, which characterizes low-lying excitations of 
a given gap structure. This N(E) gives rise to a specific 
power law temperature dependence [2j. For example, the 
line (point) node gap yields a C(T)/T ~ T(T 2 ) behavior 
in specific heat, n(T) ~ T 2 (T 3 ) for thermal conduc- 
tivity and Tj -1 ~ T 3 (T 5 ) in nuclear relaxation time at 
lower T region. This comes from the fact that the den- 
sity of states is described by a specific functional form; 
N(E) oc \E\ for line node and N(E) oc E 2 for point node, 
through which a simple power counting rule yields spe- 
cific power law indices in various quantities. Therefore it 
is decisive to precisely understand the DOS form N(E) 
in order to identify the gap structure. We have attained 
a lot of information of the pairing symmetry in various 
superconductors by this way@. 

We have been realizing, however, that there is no de- 
tailed study to investigate "the general rule" for describ- 



ing N(E) in the mixed state of type II superconductor 
for various gap structures Q. Namely although the above 
thermodynamic measurements are often performed un- 
der an applied field, which is particularly true for T\ in 
NMR, one needs to compute the precise functional form 
of N(E) in the vortex state for various gap structures. 

In the previous theoretical-work we studied only the 
full gap case and analyzed C(T) at low T 4]. However, 
this behavior should be examined with including the con- 
tribution of the gap anisotropy. Usually, experimental 
data taken under a field are often analyzed by simply 
extending the zero field results discussed above to a fi- 
nite field (H > 0), keeping the same power law of N(E) 
with a finite residual DOS Nq(H). Namely, it has been 
implicitly postulated that N(E) = N (H) + a\E\~< with 
7 = (U-shape) in the full gap, 7 = 1 in the line node 
and 7 = 2 in the point node cases. However, these forms 
of DOS spectrum should be exactly evaluated by micro- 
scopic calculation without ambiguous assumption. The 
functional form of the DOS is important when we discuss 
the physical quantities in the mixed state. 

Here we calculate the density of states N{E) aver- 
aged over a unit cell of a vortex lattice in type II su- 
perconductors with full gap, point and line node gap 
structures. Our microscopic computation takes into ac- 
count exactly quasi-particle contributions due to vor- 
tices, which is not captured by Ginzburg-Landau and 
London theories. The results are of the general form 
N(E) = N (H) + a\E\ + 0(E 2 ) for all cases, i.e. 7=1, 
showing that this singular V-shape DOS is universal in- 
dependent of the underlying original gap structure. 

We measured the DOS N(E) for two superconductors; 
NbSe 2 5] and YNi 2 B 2 C0 by performing scanning tun- 
neling spectroscopy (STS) at low temperatures. We in- 
tegrate the measured local DOS (LDOS) over a certain 
area around a vortex core. The former material is known 
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to posses an anisotropic gap without node while the lat- 
ter is speculated to be a point-like nodal gap. Thus these 
materials provide an excellent testing ground to check 
the theoretical prediction. Indeed both give rise to a V- 
shape DOS in common. This study may be useful in 
understanding vortices recently produced in ultra-cold 
Fermionic superfluids in 6 Li atomic gases)?). 

Theory: The quasiclassical theory is valid in the case 
of ^ 1, which is satisfied in almost all type II super- 
conductors, fcp is Fermi-wave number and £ is BCS co- 
herence length which is our units of the length scale. We 
introduce the pair potential A(r), the vector potential 
A(r) and the quasiclassical Green's functions g(itu n , r, k), 
r, k) and f^(\w ni r, k), where r is center of mass co- 
ordinate of the Cooper pair and the direction of the mo- 
mentum k = k/ |k|. The Eilenberger equation is given 
by 



• ( ? + ^A(r) ] \f = A(r,k) 5 , 
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A(r) / f - A*(r,k) 5 , (1) 



where g — [1 — f^f} 1 ^ 2 , Reg > 0, vf_ = vph is Fermi 
velocity, and 0o is a flux quantum 0, The applied 
field H is along z-direction. With the symmetric gauge, 
the vector potential is written as A(r) = x r + a(r), 
and an internal field h(r) is given by h(r) ~ V x a(r). 
By the numerical calculation Eq. Q is sclf-consistently 
solved with the assumption Vb(k, k ) = Vby(k)y>(k ) and 
A(r, k) = A(r)tp(k), considering the self-consistent con- 
ditions for A(r) and a(r); 

A(r) = N 2nT £ V (<p(k)f(ioj n , r, k)) t , (2) 



w„>0 



j(r) 



2nT £ i<k ff (iw n ,r,k)) t) (3) 
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where j(r) = V x V x a(r), No is density of states at 
the Fermi level in the normal state. (• • ■)(• means the 
average over the direction of k. The cut-off energy is set 
as uj c = 20T C . A is uniform gap at T — 0. We use the 
high k for the Ginzburg-Landau parameter. The angle 
dependence of the pair potential </j(k) specifies various 
gap structures with isotropic, line node or point node 
gaps. 

The LDOS at an energy E is given by 



N(E, r) = iVo(Re ff (iw„ -> E + it], r, k)) £ 



(4) 



where g is calculated by Eq. with \oj n — * E + ir/. We 
set 77 = O.OIAq. The total DOS N(E) is the spatial av- 
erage of the LDOS, i.e. N(E) = (N(E,r)) r . The self- 
consistent calculation is performed within the vortex lat- 
tice unit cell, which is divided into 81 x 81 mesh points. 
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(a) isotropic gap 

H/H c2 =0, 0.26, 0.35, 
0.44, 0.53, 1 

(b) line-node gap 
H/H c2 =0, 0.17, 0.25, 

0.33, 0.42, 1 

(c) point-node gap 
H/H c2 =Q, 0.017, 0.043, 

0.086, 0.17, 1 




-2.0 -1.0 1.0 2.0 
E 



FIG. 1: (Color online) Averaged density of states N(E) of (a) 
isotropic gap, (b) line-node gap and (c) point-node gap. H 
increases from bottom to top at E = 0. 



We assume that vortices form a triangular lattice. The 
cylindrical Fermi-surface is chosen for the s-wave and the 
line-node case. And we use the spherical Fermi-surface 
for the point-node case. 

We display the DOS N(E) for the s-wave case with 
the isotropic gap in Fig. 1(a). It is seen from this that 
the DOS has a V-shape form with a singularity at E = 0, 
namely, 



N(E)=N (H) + a \-p- + O(E 2 ) 



(5) 



where Nq(H) is the energy independent term, corre- 
sponding to the zero-energy DOS which is only a func- 
tion of the field H. This yields the Sommerfeld T- 
linear coefficient 70 (if) in the specific heat as 70 (if) = 
lim T ^o C(T)/T = (2vr 2 /3)iVo (if). We confirm that the 
sharp V-shape DOS is easily smeared by thermal broad- 
ening effect. Impurities also affect smearing of the sharp 
edge of the V-shape DOS which becomes round by fi- 
nite mean free path effect due to impurities. The linear 
slope a(H) in N(E) is field dependent; For the s-wave in 
Fig. 1(a), a(H) is small at lower fields. As if increases 
a(ff) takes a maximum around H/H C 2 ~ 0.3 and then 
becomes smaller again towards ff C 2, continuing to a flat 
DOS in the normal state. This non-monotonic behavior 
of the slope is the same also for the point node case, but 
is contrasted with that in the line node case, which are 
discussed below. 

In Fig. 1(b) we exhibit the results for the line node 
case described by the pair function <p(9,(f>) = \/2cos2#. 
At if = 0, N(E) = a\E\/A without a constant term in 
the low energy range as expected. In the vortex state this 
expression is valid except that we must add a field depen- 
dent E = value, that is, N(E) = N (H) + a\E\/A . 
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(a) Theory 




(b) Experiment 

A 



Bias Voltage 







1.0 



FIG. 2: (Color online) (a) Theoretical calculation of LDOS 
N(E, r), showing very similar sub gap peak structures to (b). 
(b) The spectral evolution N(E, r) in NbSe 2 at 0.025 T along 
109 nm length that intersects a vortex as measured by STM- 
STS. 



At first sight this result seems "apparent" because the 
original zero field DOS has a V-shape and it is simply 
shifted up-wards. However, this is not the case because 
the V-shape DOS is obtained even when the original gap 
structure is not a V-shape as already shown in the s-wave 
case. The field dependence of a(H) in the line node case 
is seen to be a monotonic decreasing function with H 
where the maximum slope occurs at H = 0. In Fig. 1(c) 
we show the V-shape DOS for the point-node case at 
H > 0. It is seen that at low fields the linear portion of 
\E\ is limited to low energies around E ~ 0. But as H in- 
creases a V-shape DOS feature manifests itself for wider 
energy ranges, so that the form in Eq. (5) is evident. 

This ubiquitous V-shape DOS can be understood in 
terms of the LDOS. We show the so-called spectral evo- 
lution for the s-wave case in Fig. 2(a) where the spectral 
weight N(E, r) in Eq. (4) is displayed in a plane of the 
distance r from the vortex core and the energy E/Aq. 
The zero-energy peak at the core site is split into two 
peaks as moving away from the core, that is, the two 
"trajectories" are given by r = (3\E\. This spectral evo- 
lution agrees with the STS observation on NbSe2 shown 
in Fig. 2(b). The total DOS is obtained by integrating 
this spectral weight spatially for each energy, namely, 

poo 

N(E) = N + drrS{r - f3\E\) = N + f3 2 \E\. (6) 
Jo 

The first term comes from the zero-energy peak at r = 0. 
This argument is applicable for the d-w&ve case where 
instead of the circular symmetry for the s-wave case 
the four fold symmetry must be taken account, which 
amounts to modifying the coefficient /3 2 , but the \E\ de- 
pendence remains unchanged. Thus it is understood that 
the V-shape DOS is universal, independent of the under- 
lying original gap structure. 

Experiment: NbSe2 is a typical anisotropic s-wave su- 
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FIG. 3: (Color online) Spectral evolution along 3 paths (black 
solid lines) for NbSe2 with respect to the vortex lattice for 
H = 0.1 T and 0.4 T. The LDOS data between the 2 hori- 
zontal white lines are used in the integration. 



perconductor with a gap ranging from 0.7 to 1.3 meV, T c 
= 7.2 K and H c2 (T = 0) = 3.2 T, whose layered struc- 
ture and Van der Waals surfaces are ideal for STM-STS 
experiments 5] . When an external field is applied to the 
c-axis, the vortices form a triangular lattice. A series of 
differential conductance spectra a(r, V) are taken along 
a line that extends radially through a vortex. Such an 
example taken at 50 mK and 0.025 T is shown in Fig. 
2(b). The overall features and numerous details of the 
LDOS agree very well with the theoretical calculations 
shown in Fig. 2(a) (also see Ref. 10). 

In order to deduce the spatially averaged DOS over a 
unit cell, the data are collected along the line paths that 
extend out in three different directions from a vortex as 
shown in Fig. 3. The normalized DOS are inferred by 
the ratio of conductance, <r(r, V)/a(r 7 V — 5 meV), and 
are shown as a color scale image. This data is weighted 
by radius and summed to emulate an integration over 
the unit cell. This is essentially a polar coordinate inte- 
gration: N(E) = J N(E,r,9)rdrd6. The resulting av- 
eraged DOS are displayed in Fig. 4(a) for two fields 
H— 0.1 T and 0.4 T along with the DOS at a zero field 
showing the s-wave gap function. A definite V-shape 
in the DOS at finite fields is evident between —0.6 and 
0.6 meV, and confirmed by the better fit attained by 
<j(V) = (To + a \V\ + (solid line) over that of a sim- 
pler parabolic form a(V) = (To + f3V 2 (dashed line) as 
shown in Fig. 4(b). Furthermore we note that the zero 
bias offset (t/<tn (cn is the normal state conductance) 
corresponding to the Sommerfeld coefficient ■y(H) Hn at 
low T nicely matches the specific heat experiments |ll| to 
within 5%. 
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FIG. 4: (a) Spatially averaged DOS N(E) over a unit cell at H=0, 0.1 T and 0.4 T for NbSe 2 . (b) Same data but over ±0.6 
mV range inside the gap. The solid lines with a + 6|V| + c| V | 2 are seen better fitting curves than the dash line with a quadratic 
form a + c'\V\ 2 . (c) DOS N(E) averaged over a unit cell at H=0, 0.07 T and 0.3 T for YNi 2 B 2 C. (d) DOS N(E) averaged 
over circular area with diameter 10 nm around a core in YN12B2C. 



Similar analyses have been performed on YN12B2C 
whose gap structure is speculated to be point-like nodes 
and <r(r, V) has been measured at 460 mK|6j. The sam- 
ples have T c = 15.6 K and H c2 (T = 0)= 8 T. When ex- 
ternal field is applied to the c-axis of this tetragonal sys- 
tem, a square vortex lattice is formed. At T — 450 mK, 
STS data are collected for various spatial points and bias 
voltages to yield the average DOS N(E). The results are 
shown in Figs. 4(c) and (d) for H = Q, 0.07 T and 0.3 
T. The conductance data sets for finite fields are seen 
to posses a V-shape dependence. Although it is difficult 
from these data to determine the precise functional form 
a(V) = (To + a\V\ v , it contains a linear term rj = 1. 

As is clear from the above argument leading to Eq. (6), 
we could see a V-shape tunneling conductance if we re- 
strict our integration region to a narrower region around 
a core. This procedure emphasizes the contribution from 
the vortex core which ultimately yields a V-shape DOS. 
In fact as shown in Fig. 4(d) which is obtained by inte- 
grating over a circular area with diameter 10 nm centered 
at a core, one can see a clear V-shape conductance curve 
for both fields. 

The ubiquitous V-shape DOS has important conse- 
quences on identifying the pairing symmetry, in particu- 
lar, its gap structure through the T-dependences of var- 
ious thermodynamic quantities as mentioned before. At 
low temperatures the T dependence is governed by the 
functional form of DOS. Using the V-shape DOS in Figs. 
1 and 4, the simple power law counting tells us the fol- 
lowings: specific heat C(T)/T — 'jo(H) + a'T, nuclear 
relaxation time Ti(T)" 1 oc T + c'T 3 , and thermal con- 
ductivity k(T) oc T + c"T 2 . Therefore, in the experiment 
under magnetic field, we can not simply assign the origin 
of the power law behavior as a line node, since the ori- 
gin may be the V-shape DOS due to the vortex states, 
not due to the line node gap structure. The If-behavior 
Ti(T) -1 ex T 3 under magnetic fields is a result of the 



spatial average. If we observe T\ outside of the vortex 
core by site selective NMR technique , we can unam- 
biguously identify the signal due to the line node. 

Simon and Lee0, and Won and Maki 1 1 derive a 
scaling law to describe T- and independences of vari- 
ous thermodynamic quantities for the line node d-wave 
systems. In particular, the latter authors explicitly find 
a scaling function based on the averaged DOS obtained 
from the Doppler shift idea|l5j whose functional form 
differs from ours. The Simon-Lee and Won-Maki scaling, 
which is quite successful for wide range of line node super- 
conductors (see for example, Ref. 16), can be improved 
and applied to other superconductors with different gap 
structures by using the obtained V-shape DOS. 

In summary, by the self-consistent quasiclassical cal- 
culation evaluating the vortex contribution exactly, 
we have demonstrated the averaged density of states 
in the mixed state is a V-shape, described by 
N(E) = N(E = 0) + a E \E\/A + 0{E 2 ). This formula 
is valid for any underlying gap structures; isotropic, point 
node or line node gaps. However the is small for the 
isotropic and point-node gap case in a low field. Two 
STM experiments on NbSe2 and YM2B2C unambigu- 
ously exhibit this behavior in their tunneling conduc- 
tance. The vortex-lattice geometry does not affect the 
V-shape and just slightly changes cce- We have dis- 
cussed several important consequences. In particular, 
when one identifies the gap symmetry by measuring ther- 
modynamic quantities, careful consideration is needed in 
the mixed state. 
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